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Abstract 

Leading-order three-nucleon forces that violate isospin symmetry are cal- 
culated in Chiral Perturbation Theory. The effect of the charge-symmetry- 
breaking three-nucleon force is investigated in the trinucleon systems using 
Faddeev calculations. We find that the contribution of this force to the 3 He - 
H binding-energy difference is given by AE^f ~ 5 keV. 
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1 Introduction 



Isospin violation [1, 2, 3] has recently been investigated in the context of Chiral 
Perturbation Theory (%PT). This powerful technique [4, 5] casts the symmetries of 
QCD into effective Lagrangians that are expressed in terms of pions and nucleons, 
which are the effective degrees of freedom of nuclear physics. These Lagrangian 
building blocks can then be combined in a systematic way to develop isospin- violating 
forces. Although most of the forces resulting from this procedure were anticipated 
and developed using phenomenological methods, new forces have also been found. 
It is the purpose of this work to complement the earlier work on isospin-violating 
two-nucleon [6, 7, 8, 9, 10, 11, 12, 13, 14] forces by calculating the leading-order 
three-nucleon isospin-violating force, which breaks charge symmetry. In addition, we 
estimate the contribution of this force to the 3 He - 3 H binding-energy difference. This 
is the first such calculation of isospin-violating three-nucleon forces 1 and completes 
the ;\TT calculations of isospin violation — both charge-independence breaking (CIB) 
and charge-symmetry breaking (CSB) — through the first three orders. 

Chiral Perturbation Theory is organized around power counting (dimensional 
analysis [15]), which allows estimates of the sizes of various mechanisms to be made 
in terms of the parameters and scales intrinsic to QCD. These scales (using Wein- 
berg power counting [4, 6, 16]) include the pion decay constant, /„- ~ 93 MeV, which 
sets the scale for pion emission or absorption, the pion mass, m n , which sets the 
scale for chiral-symmetry breaking, the typical nucleon momentum, Q ~ (which 
also determines the inverse correlation length in nuclei), and the characteristic QCD 
bound-state scale, A ~ m p , which is appropriate for heavy mesons, nucleon reso- 
nances, etc. The latter states are frozen out and do not explicitly appear in x?T, 
although their effect is present in the counter terms of the effective interactions. The 
resulting field theory is a power series in Q/A, and the number of implicit powers of 
1/A (e.g., n) can be used to label individual terms in the Lagrangian (viz., £(")). In 
this way higher powers denote smaller terms, and this is a critical part of the orga- 
nizing principle of xPT. We note that power-counting estimates of sizes are typically 
within a factor of 2-3 of the actual sizes. 

The nucleon-nucleon correlation-length scale is not relevant for one-body opera- 
tors, occurs once for two-body operators, twice for three-body operators, etc. Thus 
it is important to incorporate this mechanism into the power counting if we wish to 
compare mechanisms with differing numbers of interacting nucleons. In addition to 

1 While this manuscript was being written, we learned of a similar investigation — using, however, 
a different choice of fields — by Epelbaum, Meifiner and Palomar [37], where there is no attempt of 
a calculation of binding energies. 
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adding the indexes "n" for each of the individual Lagrangians (see above) that are 
used in a given calculation, we must add 2 for each loop and N-l for an iV-body Feyn- 
man diagram in order to determine the effective order A [16]. Thus three-nucleon 
forces with an order determined by an index n — 2 should be comparable in size in 
nuclei to one-body operators (such as the kinetic energy) corresponding to an index 
n, or a two-body force corresponding to an index n — 1. (Notice that this accounting 
of the relative sizes of few-body forces is different from Weinberg's [4, 6] by one order, 
which reflects the fact that we are counting the nucleon mass Mjv as A, rather than 
A 2 /Q.) This is the underlying reason why iV-body forces in nuclei get systematically 
smaller as N increases (and this makes nuclear physics tractable). 

Isospin violation in nuclei arises from three distinct mechanisms. The first is the 
up-down quark-mass difference, which dominates and makes the neutron heavier than 
the proton. The second mechanism is hard electromagnetic (EM) interactions at the 
quark level, which tries to make the proton heavier than the neutron. This is also the 
mechanism that produces most of the pion-mass difference. The final mechanisms are 
the soft-photon interactions (such as the Coulomb interaction between protons) that 
dominate isospin violation in nuclei. 

Direct comparison [10, 11] of the sizes of the EM and quark-mass terms demon- 
strates that the EM terms (which contain a factor of a, the fine-structure constant) 
are roughly the same size as quark-mass terms that are formally three orders smaller 
in the power counting. We adjust our power counting accordingly and adopt the 
convenient mnemonic of adding 3 to the order of the EM-induced isospin-violating 
Lagrangian when comparing sizes with quark-mass-induced mechanisms. Henceforth 
our power counting for any EM-induced interactions will contain this additional factor 
of 3. 

Our prior work on isospin-violating nucleon-nucleon forces (both CSB and CIB) 
in the context of Chiral Perturbation Theory with the A-isobar integrated out [6, 7, 
8, 9, 10, 11] identified the leading mechanisms for isospin violation in nuclei. In the 
next Section we discuss their impact on the 3 He - 3 H binding-energy difference. In the 
following Section we calculate the leading isospin-violating three-nucleon force and 
evaluate its contribution to 3 He - 3 H binding-energy difference. We then conclude. 



2 Various CSB Mechanisms 

We have shown [6, 7, 10, 11] that the following ten mechanisms are expected to 
contribute dominantly to CSB in nuclei: 
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• The mass difference of the proton and neutron, <5M N = m p — m„ < 

• The CSB nuclear kinetic energy 

• The Coulomb potential between protons 

• The Breit-interaction ((v/c) 2 ) corrections to the Coulomb potential 

• The CSB one-pion-exchange potential (OPEP) 

• The CSB short-range nuclear potential 

• The CSB two-pion-exchange potential incorporating the nucleon-mass difference 

• The Class IV CSB interactions (anti-symmetric in isospin coordinates and with 
a spin-orbit-type spin-space dependence) 

• An OPEP with Class IV isospin structure that vanishes in the two-nucleon 
center-of-mass, but not in a three-nucleon system 

• A two-pion-exchange three-nucleon force proportional to the quark-mass-difference 
contribution, 5M^ m , to the nucleon-mass difference 

We will briefly discuss each of them in turn in the context of the 3 He - 3 H binding- 
energy difference. 

The first four mechanisms are fairly well-known. They include the two largest 
mechanisms, and we start our discussion with them. 

The mass (rest-energy) difference of the nucleons, <5Mn, contributes to the xPT 
Lagrangian at order n — 1. From a nuclear-physics perspective it makes an unin- 
teresting contribution to the mass difference of 3 He and 3 H and is conventionally 
removed, leaving only a binding-energy difference. Although 3 H is heavier than 3 He, 
this removal leads to 3 He being less bound than 3 H by 764 keV, which is the target 
for all CSB calculations in the three-nucleon systems. The nucleon-mass difference 
nevertheless plays a non-trivial role in intermediate states where two protons are con- 
verted to two neutrons (or vice versa) by exchanging pions. That effect was recently 
treated in a systematic fashion [11] by removing the <5Mn mass-difference term from 
the ;\TT Lagrangian. This removes <5Mn from asymptotic states and nuclear energies, 
but its effect in intermediate states is compensated by the addition of new terms in 
the Lagrangian that must be incorporated in any calculations. The resulting scheme 
is much simpler to use than older techniques, and we will use it below. 
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The kinetic-energy difference between two protons and two neutrons caused by 
their different masses corresponds to n = 3 in power counting. In the trinucleon 
systems this mechanism leads to a robust 14 keV contribution [17, 18, 19, 20, 21, 22] 
to the binding-energy difference of 3 He and 3 H. 

The Coulomb potential between two protons is the dominant CSB interaction 
in nuclei. According to the way we bookkeep EM interactions this is an effect one 
order down compared to the leading, isospin-conserving nucleon-nucleon force, so it 
is effectively an n = 1, or A = 2, term. This contribution has a nominal size in terms 
of scales given by Eq ~ aQ ~ 1 MeV, where a is the fine-structure constant. In 
the trinucleon systems it has been well studied over several decades and leads to a 
robust and dominant 648 keV contribution to the 764 keV trinucleon binding-energy 
difference [23]. 

Small EM contributions of relativistic order contained in the Breit interaction 
(viz., the interaction between nucleon magnetic moments and between the currents 
associated with moving protons) plus a smaller vacuum-polarization force appear 
two orders down (n = 3 or A = 4 in our power counting). In terms of scales the 
relativistic contributions behave like E-q ~ aQ 3 /A 2 ~ 25 keV. Indeed, they lead to a 
robust [17, 18, 19, 20, 21] 28 keV. 

The effect of CSB on two-nucleon potentials is subsumed by the next four mech- 
anisms on the list above, of which three are Class III and one is Class IV. We will 
discuss them separately. 

Charge-symmetry breaking in the pion-nucleon coupling constants can lead to a 
CSB OPEP that has nominal size n = 2, which corresponds to A = 3. Only an 
upper limit of size 50 keV (with unknown sign) constrains this mechanism [10]. A 
conventional short-range interaction of undetermined strength corresponding to size 
A = 3 (and a nominal size of roughly 50 keV) is also present [7]. The last of the 
three Class III mechanisms is the recently calculated two-pion-exchange potential that 
incorporates various aspects of the nucleon-mass difference. It is of nominal order n = 
3 or A = 4 [10]. Each of these three mechanisms contributes to the difference between 
the pp force (with the EM interaction removed) and the nn force. At present the 
only experimental information on this difference is contained in the scattering-length 
difference. The resulting a nn —a pp scattering-length difference [2, 3, 7] of —1.5(5) fm is 
then attributed to CSB in the three forces discussed above, which cannot be further 
disentangled at the present time. (Of course, in principle these mechanisms could 
be separated thanks to their different ranges, provided that the nucleon-nucleon data 
is accurate enough.) This difference then produces a contribution to the 3 He - 3 H 
binding-energy difference of approximately 65(22) keV, a number that also appears 
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to be robust [20, 21, 24]. 



The Class IV two-nucleon CSB OPEP [2, 3, 11] has a nominal n = 2 or A = 3 size, 
which is suppressed by nearly an order of magnitude by nature's fine tuning of SM-^ to 
its physical value [11]. This type of force is further suppressed in the trinucleon bound 
states because S-wave components of the wave function do not contribute to a spin- 
orbit force. Although this force and a short-range force of order n = 4 are formally 
part of the the CSB two-nucleon potential, they make a negligible contribution to the 
trinucleon binding-energy difference. 

The last two mechanisms are three-body effects. 

The recent study [11] of Class IV CSB forces found a peculiar two-body force 
that vanishes in the center-of-mass of two nucleons, but does not vanish in a system 
of more than two nucleons. Although nominally of order n = 2 (or A = 3), this 
force should be much smaller than that for three reasons. The first reason is that 
this two-body force is constrained by kinematics to vanish in the center-of-mass of 
those two nucleons. In addition, this force is proportional to <5Mn, which results from 
the cancellation of the separate quark-mass and EM contributions and has been fine- 
tuned by nature to a rather small value (a factor of 5 smaller than the nominal value 
of the power-counting estimate for the quark-mass part SM^ m of that mass difference, 
viz., —7 MeV). The final suppression is caused by approximate SU(4) symmetry in 
the few-nucleon systems. The spin-isospin dependence of the force is antisymmetric 
under the interchange of those coordinates for the two nucleons, caused by a (rj x Tj) z 
type of isospin dependence. The dominant component of the trinucleon wave function 
(~ 90%) is the S-state (an SU(4) classification), which is completely antisymmetric 
under that interchange. These symmetry considerations cause the diagonal S-state 
matrix element of the force to vanish. The net result of these suppressions is that 
this force should be much smaller than its nominal order indicates (i.e., A = 3) and 
is therefore very unlikely to be significant. 

The remaining force is a three-nucleon force of nominal order n = 1 or A = 3 that 
originates in the chiral-symmetry-breaking properties of the quark-mass difference. 
It has never before been calculated, and we turn to it in the next Section. 



3 CSB Three-Nucleon Forces 

In this section we examine isospin-violating three-nucleon forces within Chiral Per- 
turbation Theory. Isospin-conserving three-nucleon forces have been derived within 
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this approach in Refs. [25, 26]. We follow the same method here. In particular, we 
ignore terms that cancel against recoil in the iteration of the two-nucleon potential. 



The field redefinition that we employed in Ref. [11] eliminated the nucleon-mass 
difference in the free Lagrangian at the cost of additional effective interactions pro- 
portional to powers of that mass difference. Only Lagrangian terms that had explicit 
time derivatives generated additional terms. Incorporating the results of that field 
redefinition through orders n = and n — 1 in the Lagrangian (including short-range 
two- and three-body terms) plus several other terms from Ref. [6] leads to the fol- 
lowing terms that contribute to isospin- violating three-nucleon forces at orders n — 1 
(CSB) and n = 2 (CIB), plus omitted terms that would contribute only to higher 
orders: 

+ (i) 

J TT 

The first of these terms breaks charge symmetry, while the remaining two break 
charge independence. We will focus here on the first term, which is the largest of all 
(n — 1). We discuss corrections at the end of this Section. 

Using the lowest-order isospin- conserving Lagrangian 

= I[7r2_ ( y 7r) 2_ m 2 7r 2] +7V t[^ _ _L T -(irxir)]N+^N^a-V(T-n)N+... , 

2 4/^ 2j n 

(2) 

a simple calculation along the lines of Ref. [26] leads to the following three-nucleon 
force corresponding to A = 3. We define the total three-nucleon force W as 

W = W 1 + W 2 + W 3 , (3) 

where the subscript refers to the number of the nucleon that emits both pions (the 
other two nucleons each absorb one of those pions), as shown in Fig. 1. The ex- 
pressions Wi are symmetric under the interchange of nucleons j and k. We then find 
that 

AM qm n 2 m 2 

wr = - g/4 * ( • *i* n\^\) a 3 • x 13 Hif 13 i) ) 

x (n • r 2 rf + Ti • r 3 r 2 z - r 2 • r 3 rl) , (4) 

where Tj is the isospin operator for nucleon i, &i is the spin operator for nucleon i, 
the vector from nucleon j to nucleon i, 

Y(x) = exp (— m n x)/(m n x), 



7 



k i j 



Figure 1: Leading isospin- violating three- nucleon force Wi, which is charge-symmetry 
breaking. A solid (dashed) line represents a nucleon (pion), and the cross the interac- 
tion generated by the quark- mass difference component of the nucleon-mass difference 
(first term in Eq. (1)). 

qa — 1-25 is the axial-vector constant, and 5M^ m is the currently unknown quark- 
mass portion of the nucleon-mass difference. 

The three- nucleon force in Eqs. (3,4) is charge-symmetry breaking. It is re- 
markable that this force appears formally at the same order as the leading isospin- 
conserving three-nucleon force [25, 26] in the theory without an explicit A-isobar. 
That is, the factors of Q and A are the same in both isospin-conserving and CSB 
three-nucleon forces — although, of course, the CSB force is down by a factor of 
£ = (m d - m u )/ (m d + m u ) ~ 1/3. 

CSB is thus relatively large among three-nucleon effects, an unusual phenomenon. 
It is agreed that three-nucleon effects (those not fixed by two-nucleon data) provide 
about 1 MeV of the three-nucleon binding energies. The leading three-nucleon forces 
contain relatively large sub-leading interactions (due to effects of the A-isobar), per- 
haps by a factor of 3 or so. Combining this with a factor of e, we could expect the CSB 
force to contribute as much as 100 keV to the three-nucleon binding-energy difference. 
Indeed, using the replacement /q = (qa , m n /2 f^) 2 /An = 0.075 (which is strictly valid 
only if the Goldberger-Treiman [28] relation is exact, that is, to lowest order), Eq. 
(4) can be written as — 25M^ m f^/g^ times a dimensionless function of coordinates 
(in units of spins, and isospins. Assuming that the matrix elements of this 

function give numbers of order 1 and that 5M^ m has its naive-dimensional-analysis 
value of —7 MeV, we arrive at 50 keV as an estimate for the size of the CSB three- 
nucleon force. This is significant, but obviously it could differ from the actual value 
by a factor of a few. The two sources of uncertainty in the size of this force are the 
values of 5M^ m and of the dimensionless function above. 

In order to better estimate the size of this CSB three-nucleon force, we have 



implemented it in our Faddeev codes. The cutoff parameter in the TM' force [26, 27] 
was adjusted slightly to produce the correct binding energy for 3 H when used in 
conjunction with the AV18 two-nucleon force [18]. In any numerical calculation it is 
necessary to regulate the Yukawa function, Y(x) in Eq. (4), and this was done in 
a way that is consistent with the Tucson-Melbourne force [27]. Perturbation theory 
was then used to calculate the binding-energy difference of 3 He and 3 H. We find 



77'CSB O w 



5M* m 

' MeV 



keV , (5) 



which is about a factor 3 smaller than our estimate. Any other set of realistic two- 
and three-nucleon forces should give similar results. 

As mentioned above, the actual value of 5M^ m is uncertain. It has been suggested 
[29, 30] that it could be extracted from pion-pro duct ion experiments [31, 32], but it is 
unclear if this can be achieved in the near future. It is likely to be smaller by a factor 
of a few than the naive estimate, so —7 MeV is to be viewed as an overestimate. Using 
—2.5 MeV for SM^ m , the contribution of our three-nucleon force is listed in Table I 
together with all the other significant contributions to the 3 He - 3 H binding-energy 
difference that we have discussed. 



The three-nucleon results are in agreement with experiment when the error bar 
associated with the strong- interaction CSB strength is taken into account. This 
conclusion is also consistent with the CSB results extracted in Ref. [33] for A = 6 — 10. 

There are, of course, other isospin-violating three-nucleon forces, but they are 
higher order in our power counting and thus should be smaller. Some are generated 
by sub- leading interactions, such as depicted in Fig. 2. The second term in Eq. 
(1) reflects the additional amount that should be added to the charged-pion mass 
(squared) in all pion propagators in isospin-conserving three-nucleon forces, such as 
the TM' force [26], which comes from the sub-leading isospin-conserving Lagrangian 



Table 1: Contributions to the 3 He - 3 H binding-energy difference in keV. The Coulomb 
interaction and associated (relativistic) Breit-interaction corrections dominate, while 
the CSB kinetic-energy difference (K.E.), the sum of the short-range two-body CSB 
force mechanisms, and the three-nucleon CSB force (calculated here for the first time) 
all make significant contributions. (In the three-nucleon force, we used 5M^ m = —2.5 
MeV for illustration.) "Theory" labels the sum of these mechanisms. 



Coulomb Breit K.E. Two-Body Three-body 


Theory Experiment 


648 28 14 65(22) 5 


760(22) 764 
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Figure 2: Sub-leading isospin- violating three-nucleon forces from sub-leading interac- 
tions. The cross represents the pion-mass difference (second term in Eq. (1)), while 
the circled cross stands for the sub-leading isospin-violating seagull (third term in 
Eq. (1)); a circle represents an interaction from the sub- leading isospin-conserving 
Lagrangian. 

£W. The third term in Eq. (1) is an isospin-violating contribution to the nucleon a- 
term, often called c\. It modifies the three-nucleon force that is generated by charged- 
pion exchanges in that interaction. Both of these modifications break, of course, 
charge independence, but not charge symmetry. These forces are transparently easy 
to implement, and we refrain from writing explicit forms. They correspond to n = 2 
or A = 4. 

At the same sub-leading order there are also soft-EM forces, where a pion and a 
photon are in the air at the same time, as in Fig. 3. While in flight between two 
nucleons, a charged pion can exchange a photon with the third nucleon. The photon 
couples either i) to the charge of the nucleon and the energy of the pion, which in 
the nuclear environment is Q 2 /Mn, or ii) through the momentum of the pion and the 
magnetic moment of the nucleon, which is a 1/M/v effect contained in the sub- leading 
Lagrangian In addition, there can be simultaneous emission of a photon and 

a charged pion by one nucleon followed by their absorption on two other nucleons. 
This can happen when the photon couples i) to the pion- nucleon vertex through 
the gauging of the axial- vector coupling (third term in Eq. (2)), and to the nucleon 
magnetic moment, or ii) to the nucleon charge, and to the pion-nucleon vertex through 
the gauging of the relativistic correction to the pion-nucleon coupling contained in the 
sub- leading Lagrangian These mechanisms are formally n = 2 in power counting 
(i.e., —1 + 3), and are suppressed by one power of Q/M N compared to the CSB 
three-nucleon force we calculated above. These three-nucleon forces therefore also 
correspond to A = 4. They break both charge independence and charge symmetry. 
Notice that they are entirely determined by gauge and Galilean invariance in terms 
of known parameters (the axial-vector coupling of the pion, the pion charge, the 
nucleon charge and magnetic moment, and the pion and nucleon masses). Effects 
from integrated-out resonances (most importantly the A-isobar) only appear one 
further order up. A subset of these EM effects has been calculated before [34, 35, 36]. 



10 



Figure 3: Sub-leading isospin-violating three-nucleon forces from pion-photon ex- 
change. A wavy line represents a photon and a circle an interaction from the gauging 
of the sub-leading isospin-conserving Lagrangian. 

We expect that these uncalculated parts of the CSB force corresponding to higher 
orders in the power counting contribute only a few keV or less, which is roughly the 
level of uncertainty in the EM corrections discussed above. 



4 Conclusion 



After discussing various charge-symmetry breaking mechanisms in nuclei, we derived 
the leading isospin-breaking three-nucleon force in Chiral Perturbation Theory with- 
out an explicit A- isobar, given by Eqs. (3, 4). This force is charge-symmetry breaking 
and appears formally at the same order as the leading isospin-conserving three-nucleon 
force. CSB could thus be a relatively large three-nucleon effect. Its strength depends 
on 5M^ m , the contribution from the quark-mass difference to the nucleon-mass dif- 
ference. We therefore can directly tie QCD to a three-nucleon effect. Unfortunately 
the actual value of 5M^ m has not yet been determined in a model-independent way 
from low-energy data, nor from lattice QCD. 

We have also, for the first time, calculated the contribution of this force to the 
3 He - 3 H binding-energy difference, given by Eq. (5). Taking SM^ m = —2.5 MeV 
for illustration, we find that 5 keV can be attributed to this force. This value is the 
same sign as the observed difference, and is somewhat smaller in magnitude than 
expected from naive dimensional analysis. As a consequence, it does not upset the 
agreement between theory and experiment when the uncertainty in two-body effects 
is accounted for. 
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